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We study theoretically, two doped Carbon nanotubes connected via a chemically active bond as 
the basic elementary gate for a recent proposal for qubit teleportation in the solid-state. We show 
that such a system provides the necessary entanglement between electron and electron-hole that is 
necessary for qubit teleportation. We simulate the transit of a defined qubit using time-dependent 
density functional theory within our structure. Finally, we critically discuss the implications of our 
design. 

PACS numbers: Valid PACS appear here 



X3 

o 



> 
in 

00 • 

o ■ 
in : 
o . 

a • 

S; 
^ . 

o . 

Q ■ 



X 



BACKGROUND 

For computers to become faster and faster, quantum 
computing must replace or supplement what we already 
haveQ. To this end, there have been many proposals 
for structures capable of executing quantum computing. 
Proposed quantum computers utilise quantum degrees of 
freedom to create quantum bits (qubits). The intrinsic 
nature of the qubit allows for the performance of tasks 
such as factoring large integers and searching for and 
simulating quantum mechanical systems. However, the 
requirements of quantum computing present many new 
and unique challenges Q. Many of the current quantum 
computing proposals rely on the optical manipulation of 
information but they also have solid-state counterparts. 

A qubit is any system that can assume either of two 
super-posed values. It is normally composed of a two 
state system where the two states are simply called |0) 
and The difference between a (classical) bit and a 
qubit being that the qubit lacks the characteristic (classi- 
cal) energy barrier separating the states hence the qubit 
is neither fully |0 or |1 but a superposition of both. The 
qubit can be written as 



qubit = ao|0) + ai|l) 



(1) 



where |0) denotes the state in which the qubit has a value 
of and |1) denotes the state in whichthe qubit has a 
value of 1. The coefficients ao and ai are complex quan- 
tities whose squared magnitudes denote the probability 
that if a measurement is performed on the qubit, it will 
be found to have a value of or 1 respectively. 

A ready example of a qubit capable system is the spin- 
up or spin-down character of the electron. To date, many 
of the solid-state schemes proposed utilise the Nuclear 
spin. Naturally the spin of an electron would be the most 
natural candidate for a qubit but visualising a system 
capable of utilising the electron in this way often involves 
optical parallels. In which case the concept of a free 
electron is employed. To date, there is no such thing. 



Whilst qubit teleportation has been demonstrated to 
be possible for optically based schemes, there is still no di- 
rectly viable occasion of qubit teleportation in the solid- 
state. The present day consensus is to move towards 
systems that can exhibit quantum computing character- 
istics and from there towards fully capable systems. In 
the solid-state case; this involves the demonstration of 
qubit teleportation. 

A design for teleportation of the charge of an electron 
(and thus the qubit) within a "Fermi-sea" has been pro- 
posed by Prof. Beenakker and his groupQ. A full sum- 
mary is available from In the following section we 
summarise their approach. In realising this scheme we 
confirm and expand upon the results of Gabor et al.Q- 
Our aim in this paper is to show how this scheme might 
be realised using current technology in the form of net- 
works of connected carbon nanotubes doped with impuri- 
ties. In realising quantum computation; the present chal- 
lenge is to move towards systems where quantum com- 
putation phenomena can be demonstrated and controlled 
in experiment. It is with this view that we present this 
paper. 



INTRODUCTION 

Electron tunnelling is a manifestation of the Heisen- 
berg Uncertainty principle and quantum mechanics. An 
electron hole is an unoccupied state within a Fermi-sea. 
It is one of only two charge carriers in the solid-statefisj 
The opposing charges of hole and electron means that 
they can only exist with a separating insulating barrier 
between them. If this barrier passes a small current and 
an electron meets a hole, then both are annihilated. The 
proposed scheme works with the use of electron and hole 
pairs separated by potential barriers, Figure^is a simple 
schematic. 

At low temperatures, all energy levels are filled with 
electrons to the maximum energy with the maximum en- 
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ergy being the Fermi-energy. This is referred to as the 
Fermi sea. This 'sea' is in equihbrium and thus carries 
no current (excitations may induce a current). The tun- 
nelling event happens when the Fermi energy on one side 
of the barrier is higher than that on the other side of the 
barrier. So the hole is raised to the same energy as the 
electron. There is a probability that the electron and the 
hole are reflected by the barrier (depending on the condi- 
tions) but each time they meet, there is also a probability 
that the electron will tunnel through the barrier and fall 
into the hole on the other side. This mechanism forms 
the basis of the teleport proposal. 

In its simplest form, the teleportation process begins 
with the entangled state (|T)e|T)fe)/'\/20. The subscripts 
e and h refer to the electron and the hole at two different 
locations respectively. The particle teleported is another 
electron in the state a|t)e + /3|i)e with (jap + = 1). 
The second electron e can tunnel into the hole h only if 
the spins match. The probability of this happening can 
be expressed as = ^laPl^P + ^l/^PKP < 1 where 
t is the tunnelling amplitude. There is no instantaneous 
transfer of information. The unpredictability of the suc- 
cess of the tunnelling attempt means that a message will 
need to be sent by classical means that the teleportation 
has happened. 

The instantaneous transfer of the state from one to 
the other electron is necessary to satisfy the no-cloning 
theorem of quantum mechanics At no point is there 
more than one copy of the state. As usual, all the limi- 
tations of teleportation still applyP,0]. We should note 
that whilst qubit teleportation has been demonstrated 
to be possible for optically based schemes, there is still 
no directly viable occasion of qubit teleportation in the 
solid-state. 

Central to the design of this 'qubit-teleporter' is the 
architecture that allows for the production (and detec- 
tion) of entangled electron-hole pairs in the Fermi sea. 
Figure „ (from the same group) shows how this might be 
achieved. This would be the building block of any solid- 
state solution for the proposed scheme. 
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THE SYSTEM 

Our aim in this paper is to highlight a system capable 
of this entanglement (and ultimately qubit teleportation) 
in the form of impurity doped carbon nanotubes con- 
nected via chemically active bonds. Previous studies 0,0 
have shown that suitably doped semi-conducting car- 
bon nanotubes form chemically active C-C bonds when 
placed in proximity to each other. These bonds can serve 
as the necessary tunnel barriers depicted in Fig ^ thus 
allowing for its realisation. 

Previous studies have also indicated that semiconduct- 
ing carbon nanotubes with chirality (8,0) are well suited 
to such a role. LDA band structure calculations predict 
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FIG. 1: Depiction of the teleportation scheme from 0. A 
voltage V is applied over a tunnel barrier produces pairs of 
entangled electron-hole pairs within the system. (eL,h_t) are 
such a pair. As in we assume spin entanglement in the state 
( ! tr ) + 1 U ) ) 7^2 where the first arrow refers to the electron spin 
and the second arrow refers to the hole spin. A second electron 
Br is in an unknown state a|t) -I- The electron e_R can 

annihilate with the hole h_L by tunnelling through the barrier 
at the centre. If it happens and is detected, then the state e^ 
collapses to the state e_R. Note that |t) annihilates with 
and II) annihilates with |j) so e^ inherits the coefficients a 
and /3 of en after its annihilation. The diagram shows a second 
entangler at the right, to perform two way teleportation. This 
essentially leads to entanglement swapping, sl and h_R become 
entangled after the annihilation of h_L and e/j 



a minimum band gap of 0.6 eV whilst Quasiparticle cor- 
rections open this up to 1.75eV 6]. Quasiparticle band 
structure calculations have shown nanotubes with chiral- 
ity (8,0) to be of an LDA determined 1.75 eV with a 
stable structure diameter of j Inm. 

The nodes in the circumferential direction within this 
structure are also quantised. With the inclusion of impu- 
rities (via doping), carbon nanotubes exhibit a reduced 
band gap and excess spin-density localisation 0, Q lo- 
calised about the impurity. The studies of have also 
indicated that the confinement of electrons and holes 
within the un-doped structure is of the order of several 
A. Within our planned structure this has the effect of 
increasing the probability of electron tunnelling and en- 
tanglement swapping. 

We simulate 2 doped semiconducting carbon nan- 
otubes connected via a chemically active bond as our 
'quantum entangler'. The tubes are of chirality (8,0) and 
doped with Nitrogen (N) . The tunnel barrier between the 
sections of nanotube is analogous to the split gate elec- 
trode in Fign The appHcation of a voltage across either 
of the tubes (across the tunnel barrier) then results in 
tunnelling of an electron on one side to an empty state 
on the other side. In our calculations we simulate 2 sec- 
tions both 2 unit cells long in the Z direction and each 
with a nitrogen impurity. 

Regard to the information contained on the tunnelling 
electron means that the barrier must be shielded from 
interference for its surrounding environment. A recent 
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FIG. 2: Schematic description of a method to produce and 
detect entangled edge channels in the quantum Hall effect 
from Q. The bold black lines mark the boundaries of a two- 
dimensional electron gas. A strong perpendicular magnetic 
field B ensures that the transport near the Fermi-level Ef 
takes place in two edge channels, extended along a pair of 
equipotentials (thin solid and dashed lines, with arrows that 
give the direction of propagation). The split gate electrode 
(shaded rectangles at the center) divides the conductor into 
two halves, coupled by tunnelling through a narrow opening 
(dashed arrow). If a voltage is applied between the two halves, 
then there is a narrow energy range eV above Ef in which 
the edge channels are predominantly filled in the left half 
(solid lines). Tunnelling events introduce filled states in the 
right-half (black dots) and empty states in the left half (open 
circles). The entanglement of these particle-hole excitations 
is detected via the violation of a Bell inequality. 

study of supcr-luminal electron tunnelling suggests that 
the information of a tunnelling electron wave-packet is 
essentially readjusted within the barrier^ |^. In effect, 
the dwell time of the electron |Q| must be kept to a min- 
imum to ensure the integrity of the wavepacket (which 
is the spin character of the electron). Furthermore, to 
ensure entanglement, the material should be a semi- 
conductor to preserve the initial state of the qubit and 
allow for the needed decoherence times. Our choice of 
the (8,0) carbon nanotube meets all these criteria. 

Our design exploits the widely reported and verified 
properties of carbon nanotubes. Within the configura- 
tion laid out above, the electron and hole are both con- 
fined to a 1 dimensional line. The only degree of free- 
dom allowed is through the tunnel barrier connecting the 
tubes. Any other attempt to escape the structure incurs 
a massive energetic penalty. In the case of carbon nan- 
otubes (with a band gap of j 2ev) the potential is still 
enough to confine the electron and hole even at room 
temperatures. In the case of a defined qubit (which is 
weakly coupled to its environment), the spatial degree 
of freedom allowed by the tunnel barrier heightens the 
tunnelling probability. 

The use of semiconducting carbon nanotubes in par- 
ticular offers many advantages. Among these are the fact 
that advances in fabrication are making the production 
of carbon nanotubes less and less expensive. There is also 
the fact that carbon nanotubes are very good structures 
for conserving the spin of the charge carrier, thus pre- 
serving the integrity of the qubit over the length scales 




FIG. 3: Classical depiction of two sections of semiconduct- 
ing Carbon nanotubes connected by chemically active bonds. 
Both sections of tube are of chirality (8,0) and 2 unit cells 
in length. The blue spheres denote impurity sites (Impurities 
being Nitrogen atoms). The excess electrons that the impu- 
rities carry are effectively localised within the structure^ 

involved. Another advantage is that one then has the 
option of layout and circuit design and all this is within 
the reach of current technology. An illustration of our 
proposed design is Fig|21 
Figure El 

METHOD 

Within our simulation cell, we define of two doped sec- 
tions of carbon nanotube with chirality (8,0) and com- 
posed of 2 primitive unit cells (each primitive unit cell 
contains 32 atoms) in the z direction. The atomic co- 
ordinates are generated using the program Tubegenj^ 
with a C-C bond length of 1.425A(as confirmed by ex- 
periments). Each section has a single nitrogen impurity 
in place of a C atom and both are in the x-plane facing 
each other in the middle as shown in Figure El There 
is an inter-tube distance of 1.6A(half the interlayer dis- 
tance in natural graphite). 

We define an excess spin density in the form of an elec- 
tron with a spin-up character. We assign this excess spin 
density to the tube section to the left in FigureEl as indi- 
cated in Figure^l Due to computational constraints, we 
only simulate 2 sections of carbon nanotube as opposed 
to infinite length tubes. However, with our approach, 
our results are easily applicable to the infinite length 
tube case. To account for our finite system we perform 
total ener gy c alculations using the Louis Bearends (LB) 
functional [ill. The LB functional is of a form designed to 
get the correct asymptotic behaviour; yielding good en- 
ergies and excitation energies in calculations of response 
properties. 
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We use improved troullier-martins pseudopotentials 
within the GGA approximation as implemented in the 
octopus code0. The pseudopotentials are distributed 
with the package and they have been extensively tried 
and tested in many different applications. We perform 
static total-energy calculations and find that when these 
localised impurities face each other they form a chemi- 
cally active bond connecting the two sections of carbon 
nanotube together (confirming the results of 

We then compute the time-dependent simulation of a 
polarised electric field across the tubes. The electric field 
is of delta strength (voltage) O.OSeV— 1 polarised -100% 
in the x direction and 0% in the y and z directions respec- 
tively. We propagate the time-dependent Kohn-Sham 
equations using the magnus evolution operator [l^. at a 
time step of O.OOOSev for a total of 2000 time steps. From 
a multipole representation of the time propagation, we 
calculate the evolution of the dipole-strength function of 
the system. In experiment, the optical-absorption cross 
section is directly accessible but in this instance, our use 
of the dipole strength function is of a qualitative nature. 
It presents a ready picture of the shifts in polarizabil- 
ity across the system as well as serving as a means of 
identification of such a system. 

Our use of the dipole strength function serves as a poor 
man's measure of the electron hole interaction. Previ- 
ous studies have shown that the optical response of 
(8,0) carbon nanotubes is dominated by excitonic terms. 
Contained within the qualitative expression of the dipole 
strength function are the excitonic effects between the 
electron, the hole and the hole created by the transition 
of the electron across the tunnel barrier. The changing 
values (and sign) of the dipole strength function with 
respect to energy is indicative of the electron-hole inter- 
action. 

To calculate the dipole-strength function (which is 
simply proportional to the absorption cross section), we 
first calculate the ground-state energy and then excite 
the system with an (afore-mentioned) delta electric field, 
EoS{t). The dipole strength can then be related to to the 
imaginary part of the dynamical polarizability by 

S{uj) = -^Lojaiio) (2) 

where h is Planck's constant, rrie the electron mass and 
the dynamical polarizability, 

a{uj) = — / drZSn{r,uj). (3) 

In the last expression, (5n(r,a;) stands for the Fourier 
transform of n(r, i) — n(r, i — 0). We present our results 
below. 



RESULTS 

From the initial static total energy calculation of the 
system, we find that the maximal force on any atom in 
the simulation cell to be less than 150(eV/A). This may 
be irrelevant given the inaccuracy of the LB functional 
in calculating forces but the figure is still within a mag- 
nitude of what is acceptable for our design. The direct 
substitutionality of the impurity atoms is confirmed; the 
maximum forces acting on the impurity atoms are less 
than 10(eV/A). We also find that at an intertube dis- 
tance of 3.2 A, (the distance between graphite layers in 
nature) the doped sections of tube form one main inter- 
tube C-C bond as opposed to the classical picture shown 
in Figure 13 A cross section of the electron localisation 
function (ELF) of the system is shown in Figure b It 
can be seen that about midpoint along the length (in 
the z-direction) of both tubes, there is a maximum elec- 
tron localisation value 0.8 (representing the bond) and 
nominal values of 0.5 on either side of the bond repre- 
senting spin-density accumulation. Analysis of the final 
Spin-Density configuration of the system shows that the 
applied electronic excitation in the form of the applied 
E-field propels this excess spin density across the bond. 
This is mapped out in the plots of the Un-paired spin- 
density (Figures 01a),(b) and (c)) at different timesteps. 

From a static representation of our system (with a de- 
fined excess qubit) we generate a picture where we have 
a qubit value (a clearly defined spin-up character to the 
spin-density and an ELF function of the chemical bond 
in Fig0Ja) . The effect of the applied field is to propel this 
excess spin density from the left side of the system across 
the bond to the right side of the system as shown in Fig- 
ure ^c) . Analysis of the dipole strength function shows 
shifts in the orientation of the dynamic polarizability of 
the system. This is in part due to the geometric shifts in 
the position of our defined qubit across the system. 

DISCUSSION 

The bond formed between the tubes is covalent and 
thus not easily broken. With an ELF value of 0.9, the 
bond is prevalent even at temperatures elevated beyond 
the ground state. We should note that the above exam- 
ple is merely illustrative of what is possible, this same 
scheme as applied to infinite length tubes in networks 
formed by chemically active bonds favours the telepor- 
tation proposal more so. In the case of infinite length 
tubes there would be even more of a localisation and iso- 
lation of the qubit wave-packet 0] ■ We note that in 
real applications, the tunnelling probability of the qubit 
would be very low. However, the system described above 
is easily scaleable to suit the design depicted in Figure 
n The optical properties of nanotubes with chirality of 
(8,0) is strongly dominated by excitonic effects. This 
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FIG. 4: The Un-paired Spin-density of the system projected 
onto the x-axis. In (a) we see a clearly defined qubit. (b) 
shows the propagation of this qubit across the tunnel barrier, 
(c) represents the position of the qubit post time propagation. 
Throughout it can be seen that the opposing and resultant 
electron-hole is delocalised across the structure. The electron 
is localised throughout. 



feature also bodes much for our proposal. Within the 
setup outhncd above the availability of excitons on either 
side of the tunnel barrier greatly increases he probability 
of a tunnelling event and hence entanglement swapping. 
However, we should also note the the effects of electron 
phonon coupling also have a part to play. 

One of the said requirements of a system capable of 
quantum computing is one which is easily scalable with 
well characterised qubits. We believe our design to be 



FIG. 5: A cross section of the electron localisation function of 
the system from the static total energy calculation projected 
onto the z-axis. Green areas indicate areas of lowest local- 
isation whilst red areas indicate the maximum. There is a 
maximum ELF value across the bond of 0.9. 



a ready example of such a system. In implementation, 
this allows for the expression of the spin 1/2 value of the 
electron as 



IT)^|i) 



(4) 



In an experimentally viable setup, we expect that any 
device would be of the form of Fig ?? An arrangement 
such a this would have a shear modulus of 80 Gpa or 
morej3- 

As concerns a traversal time or a tunnelling time, a 
true definition is (not or might not ever) be possible 
(an un-controversial definition still escapes the litera- 
ture.) However, within the context of this work, the 
'view-point' of the qubit is what we are interested in and 
what we 'impose' on the simulation. True electron tun- 
nelling (without distortion) is of itself a quasistatic pro- 
cess in which the output density (from the tunnel barrier) 
adiabatically follows the input with a small delay due to 
energy storage within the barrier This time-delay is 
(probably) the most representative of what we might de- 
fine as an operation time for the device, or a dwell time 
for the information carrier. In this simulation, this time 
was of the order of 2eV~^ (at an intertube distance of 
1.6A). 

In general, the scheme as approximated in this model 
can only be a guide as to the decoherence times one might 
expect of the defined qubit. As a rough guide, we em- 
ploy the equation of Burkard and Loss; in which they 
posit that the spin coherence time may be measured as 
a function of the Electron Spin Resonance (ESR). The 
resulting stationary current would exhibit a resonance 
whose line width is determined by this time. The ulti- 
mate feasibility of such a configuration must be extracted 
from experiment; our results clearly support such a con- 
figuration. 
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FIG. 6: The dipole strength function of the system over 
a time-interval representative of the qubit transition time. 
A Comparison with Energy indicates multiple peaks and 
troughs. The first of which we define as the deterministic 
point of tunnelling. This generated optical signature may 
serve as an aid to the detection of the event within a bundle 
of nanotubes. 

CONCLUSIONS 

In conclusion, We have performed total-energy pseu- 
dopotential calculations on two sections of doped car- 
bon nanotube using the plane-wave pseudopotential code 
octopus^3- We have shown that when the impurities 
face each other, a bond forms between the two sections of 
carbon nanotube (confirming the results of We have 
shown that this chemically active bond may be used as 
a tunnel barrier in the realisation of the solid-state en- 
tanglement scheme of Finally, we have shown that 
such a system allows for electron hole entanglement that 



makes the solid state teleportation scheme depicted in 
work. Extending this, we have identified an optical 
profile in the form of a qualitative dipole response which 
may be used to identify and isolate such a device within 
a bundle of tubes. 
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